BABAR-CONF-02/07 
SLAC-PUB-9229 
May, 2002 



Measurements of Branching Fractions and CP- Violating 
Asymmetries in B° — > 7r + 7r _ , Ktt. K + K~ Decays 



The BABAR Collaboration 
February 7, 2008 

Abstract 



We present updated measurements of branching fractions and CP-violating asymmetries for neutral 
B meson decays to two-body final states of charged pions and kaons. The results are obtained from 
a data sample of about 60 million T(45) — > BB decays collected between 1999 and 2001 by the 
BABAR detector at the PEP-II asymmetric-energy B Factory at SLAC. The sample contains 124]ti5 
7T7T, 403 ± 24 Ktt, and 0.6^74 KK candidates, from which we measure the following quantities: 
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= -0.02 ±0.29 ±0.07 [- 
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where the errors are statistical and systematic, respectively, and the asymmetry limits correspond 
to the 90% confidence level. These results are preliminary. 
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Recent measurements of the CP-violating asymmetry parameter sin2/3 by the BABAR [T[ and 
Belle || collaborations established CP violation in the B° system. These measurements, as well 
as an updated result by BABAR Q reported at this conference, are consistent with the Standard 
Model expectation based on measurements and theoretical estimates of the elements of the Cabibbo- 
Kobayashi-Maskawa [Q] (CKM) quark-mixing matrix. 

The study of B decays to charmless hadronic two-body final states will yield important infor- 
mation about the remaining angles (a and 7) of the Unitarity Triangle. In the Standard Model, 
the time-dependent CP-violating asymmetry in the decay B° — > tt + tt~ is related to the angle a, 
and ratios of branching fractions for various tttt and Kir decay modes are sensitive to the angle 7. 
In this paper, we update our previous measurements of branching fractions ||] and CP-violating 
asymmetries Q in B° — > 7r + 7r~, K + tt^ , and K + K~ decays^ using a sample of 60 million BB pairs. 

We reconstruct a sample of B mesons (B rec ) decaying to the h + h'~ final state, where h and 
b! refer to tt or K, and examine the remaining charged particles in each event to "tag" the flavor 
of the other B meson (B tag )- The decay rate distribution / + (/_) when h + h'~ = 7r + ir~ and 
?tag 



B tag = B° (B°) is given by 



e -|At|/r 

f±{At) = — [1 ± sm(Am d At) t C nn cos{Am d At)], (1) 

4r 

where r is the mean B° lifetime, Am,; is the eigenstate mass difference, and At = t rec — tt ag is the 
time between the B rec and B tag decays. The CP-violating parameters S nn and are defined as 

_ 21m\ A n 1 - |A| 2 
If the decay proceeds purely through the b — > uW~ tree process, then A is given by 

X{B ^ vr+7r - ) = ( Ym^l] I YMA . ( 3) 
1 ; \v tb v* d ) \v ud v:j { ' 

In this case C^ = and S W7r = sin2a, where a = arg [— ^dV^/V^V^]. In general, the b — > dg 
penguin amplitude modifies both the magnitude and phase of A, so that C nn 7^ and S n7T = 
\fl — C% n sin2a e ff, where a e g depends on the magnitudes and relative strong and weak phases of 
the tree and penguin amplitudes. Several approaches have been proposed to obtain information on 
a in the presence of penguins @. 

The data sample used in this analysis consists of 55.6 fb , corresponding to (60.2 ±0.7) million 
BB pairs, collected on the T(4S) resonance with the BABAR detector at the SLAC PEP-II storage 
ring between October 1999 and December 2001. A detailed description of the detector is presented 
in Ref. ||]. Charged particle (track) momenta are measured in a tracking system consisting of a 
5-layer double-sided silicon vertex tracker (SVT) and a 40-layer drift chamber (DCH) filled with a 
gas mixture of helium and isobutane. The SVT and DCH operate within a 1.5 T superconducting 
solenoidal magnet. Photons are detected in an electromagnetic calorimeter (EMC) consisting of 
6580 CsI(Tl) crystals arranged in barrel and forward endcap subdetectors. The flux return for the 
solenoid is composed of multiple layers of iron and resistive plate chambers for the identification 
of muons and long-lived neutral hadrons. Tracks from the B rec decay are identified as pions or 
kaons by the Cherenkov angle 9 C measured with a detector of internally reflected Cherenkov light 
(DIRC). 

1 Unless explicitly stated, charge conjugate decay modes are assumed throughout this paper. 
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Event selection is identical to that described in Ref . || . Candidate B rec decays are reconstructed 
from pairs of oppositely-charged tracks forming a good quality vertex, where the B Tec four-vector 
is calculated assuming the pion mass for both tracks. We require each track to have an associated 
6 C measurement with a minimum of six Cherenkov photons above background, where the average 
is approximately 30 for both pions and kaons. Protons are rejected based on 9 C and electrons are 
rejected based on dE/dx measurements in the tracking system, shower shape in the EMC, and 
the ratio of shower energy and track momentum. Background from the reaction e + e~ — ► qq (q = 
u,d,s,c) is suppressed by removing jet-like events from the sample: we define the center-of-mass 
(cm.) angle 9s between the sphericity axes of the B candidate and the remaining tracks and 
photons in the event, and require |cos#s| < 0.8, which removes 83% of the background. The total 
efficiency on signal events for all of the above selection is approximately 38%. 

Signal decays are identified kinematically using two variables. We define a beam-energy substi- 



tuted mass m-Es = J E^ — p 2 ^, where the B candidate energy is defined as E^ = (s/2 + Pi • ps)/^, 
\/~s and Ei are the total energies of the e + e~ system in the cm. and laboratory frames, respectively, 
and Pi and p# are the momentum vectors in the laboratory frame of the e + e~ system and the 
-B rec candidate, respectively. Signal events are Gaussian distributed in mES with a mean near the 
B mass and a resolution of 2.6MeV/c 2 , dominated by the beam energy spread. The background 
shape is parameterized by a threshold function Q with a fixed endpoint given by the average beam 
energy. 

We define a second kinematic variable AE as the difference between the energy of the B rec 
candidate in the cm. frame and -y/s/2. Signal tttt decays are Gaussian distributed with a mean 
value near zero. For decays with one (two) kaons, the distribution is shifted relative to tttt on 
average by —45 MeV (—91 MeV), respectively, where the exact separation depends on the laboratory 
momentum of the kaon(s). The resolution on AE is approximately 26 MeV and is validated in large 
samples of fully reconstructed B decays. The background is parameterized by a quadratic function. 

Candidate h + h'~ pairs selected in the region 5.2 < ot-es < 5.3GeV/c 2 and \ AE\ < 0.15 GeV 
are used to extract yields and CP-violating asymmetries with an unbinned maximum likelihood fit. 
The total number of events in the fit region satisfying all of the above criteria is 17585. 

To determine the flavor of the -Btag meson we use the same B-tagging algorithm used in the 
BABAR sin2/3 analysis [1C]. The algorithm relies on the correlation between the flavor of the b 
quark and the charge of the remaining tracks in the event after removal of the B Tec candidate. We 
define five mutually exclusive tagging categories: Lepton, Kaon, NT1, NT2, and Untagged. Lepton 
tags rely on primary electrons and muons from semileptonic B decays, while Kaon tags exploit 
the correlation in the process b — > c — > s between the net kaon charge and the charge of the b 
quark. The NT1 (more certain tags) and NT2 (less certain tags) categories are derived from a neural 
network that is sensitive to charge correlations between the parent B and unidentified leptons and 
kaons, soft pions, or the charge and momentum of the track with the highest cm. momentum. The 
addition of Untagged events improves the signal yield estimates and provides a larger sample for 
determining background shape parameters directly in the maximum likelihood fit. 

The quality of tagging is expressed in terms of the effective efficiency Q = J2c € c^c^ where 
e c is the fraction of events tagged in category c and the dilution D c = 1 — 2w c is related to the 
mistag fraction w c . Table |I| summarizes the tagging performance in a data sample -Bfl av of fully 
reconstructed neutral B decays into D^~h + (h + = ir + , p + , af) and J/i/jK*° (K*° — > K + ir~) flavor 
eigenstates. We use the same tagging efficiencies and dilutions for signal nir, Kir, and KK decays. 
Separate background efficiencies for each species are determined simultaneously with S n7T and C n7V 
in the maximum likelihood fit. 



9 



Table 1: Tagging efficiency e, average dilution D = 1/2 (D B o + Dg ), dilution difference AD = 
D B o — Dg , and effective tagging efficiency Q for signal events in each tagging category. The values 
are measured in the -Bfl av sample. 



Category 




D{%) 


AD (%) 


Q(%) 


Lepton 


11.1 ±0.2 


82.8 ± 1.8 


-1.2 ±3.0 


7.6 ±0.4 


Kaon 


34.7 ±0.4 


63.8 ± 1.4 


1.8 ±2.1 


14.1 ±0.6 


NT1 


7.6 ±0.2 


56.0 ±3.0 


-2.7 ±4.7 


2.4 ±0.3 


NT2 


14.0 ±0.3 


25.4 ±2.6 


9.4 ±3.8 


0.9 ±0.2 


Untagged 


32.6 ± 0.5 








Total Q 








25.0 ±0.8 



The time difference At is obtained from the measured distance between the z positions of the 
B Tec and -B tag decay vertices and the known boost of the e + e~ system. The z position of the -B tag 
vertex is determined with an iterative procedure that removes tracks with a large contribution to 
the total x 2 - An additional constraint is constructed from the three-momentum and vertex position 
of the B rec candidate, and the average e + e~ interaction point and boost. For 99.5% of candidates 
with a reconstructed vertex the r.m.s. Az resolution is 180 /im (1.1 ps). We require \ At\ < 20 ps and 
f At < 2.5 ps, where a At is the error on At. The resolution function for signal candidates is a sum of 
three Gaussians, identical to the one described in Ref. 0], with parameters determined from a fit to 
the -Bflav sample (including events in all five tagging categories). The background At distribution 
is parameterized as the sum of an exponential convolved with a Gaussian, and two additional 
Gaussians to account for tails. A common parameterization is used for all tagging categories, and 
the parameters are determined simultaneously with the CP parameters in the maximum likelihood 
fit. We find that 86% of background events are described by an effective lifetime of about 0.6 ps, 
while tails are described by 12 (2)% of events with a resolution of approximately 2 (8) ps. 

Discrimination of signal from light-quark background is enhanced by the use of a Fisher discrim- 
inant T IJ. The discriminating variables are constructed from the scalar sum of the cm. momenta 
of all tracks and photons (excluding tracks from the B iec candidate) entering nine two-sided 10- 
degree concentric cones centered on the thrust axis of the B rec candidate. The distribution of T for 
signal events is parameterized as a single Gaussian, with parameters determined from Monte Carlo 
simulated decays and validated with B~ — > D°ir~ decays reconstructed in data. The background 
shape is parameterized as the sum of two Gaussians, with parameters determined directly in the 
maximum likelihood fit. 

Identification of h + h'~ tracks as pions or kaons is accomplished with the Cherenkov angle 
measurement from the DIRC. We construct Gaussian probability density functions (PDFs) from the 
difference between measured and expected values of 9 C for the pion or kaon hypothesis, normalized 
by the resolution. The DIRC performance is parameterized using a sample of D* + — > D°ir + , 
D° —* K~Tt + decays, reconstructed in data. The typical separation between pions and kaons varies 
from 8a at 2GeV/c to 2.5<7 at 4GeV/c, where a is the average resolution on 9 C (Fig. [j]). 

We use an unbinned extended maximum likelihood fit to extract yields and CP parameters from 
the -B rec sample. The likelihood for candidate j tagged in category c is obtained by summing the 
product of event yield rtj, tagging efficiency ej )C , and probability 7\ c over the eight possible signal 
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Figure 1: Variation of the separation between the kaon and pion Cherenkov angles with momentum, 
as obtained from a control sample of D* + — > D°-k + , D° — > K~ir + decays reconstructed in data. 



and background hypotheses i (referring to tttt, K + ir , K ir + , and KK decays), 



exp 



E 



n 



cPi. 



(4) 



For the K^ir^ components, the yield is parameterized as = Nk-k (1 ± Ak-k) /2, where A^^^ = 
^k--k++^k+tt- and^KTr = (A r x-7r+ _ ^ r ftr+7r-)/(^ r x-7r++^ r K+7r-)- Tne probabilities 7\ c are eval- 
uated as the product of PDFs for each of the independent variables Xj = {m^s, AE, J 7 , 6+,9~ , At}, 
where 0+ and 6~ are the Cherenkov angles for the positively and negatively charged tracks. We use 
the same PDF parameters for and 9~. The total likelihood C is the product of likelihoods for 
each tagging category and the free parameters are determined by minimizing the quantity — ln£. 

In order to minimize systematic error on the branching fraction measurements, we perform 
an initial fit without tagging or At information. A total of 16 parameters are varied in the fit, 
including signal and background yields (6 parameters) and asymmetries (2), and parameters for 
the background shapes in rngs (1), A£7 (2), and J- (5). Table Q summarizes results for signal 
yields, total efficiencies, branching fractions, and Ak-k- The upper limit on the signal yield for 
B° — ► K + K~ is given by the value of n° for which J Q n C maiX dn/ J °° C m8bX dn = 0.90, where >C max 
is the likelihood as a function of n, maximized with respect to the remaining fit parameters. The 
branching fraction upper limit is calculated by increasing the signal yield upper limit and reducing 
the efficiency by their respective systematic errors. The dominant systematic error on the branching 
fraction measurements is due to uncertainty in the shape of the 9 C PDF, while the dominant error 
on Ak-k is due to possible charge bias in track and 6 C reconstruction. All measurements reported 
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Table 2: Summary of results for total detection efficiencies (Eff), fitted signal yields N$, measured 
branching fractions B, and the Kir charge asymmetry Ak-k- The sample corresponds to (60.2 ±0.7) 
million BB pairs produced, where equal branching fractions for T(AS) — > B°B° and B + B~ are 
assumed. The statistical and systematic errors on Akk are added in quadrature when calculating 
the 90% confidence level (C.L.). 



Mode 


Eff (%) 


N s 


H(10- e ) 


Ak-k 


Ak-k 90% C.L. 


7T + 7T~ 


38.5 ±0.7 




5.4 ±0.7 ±0.4 






K + ir- 


37.6 ±0.7 


403 ± 24 ± 15 


17.8 ±1.1 ±0.8 


-0.05 ±0.06 ±0.01 


[-0.14, +0.05] 


K+K- 


36.7 ±0.7 


0.61?;° (< 15.6) 


< 1.1 (90% C.L.) 







in Table ^ are consistent with our previous results reported in Ref. ||. 

Figure [2] shows distributions of tbes an d AE after a cut on likelihood ratios. We define 
7£ s ig = J2 S n sP s l J2i n{Pi and IZ^ = n^Vk/Yls n s'Ps, where J2s (Xi) indicates a sum over signal (all) 
hypotheses, and indicates the probability for signal hypothesis k. The probabilities include the 
PDFs for 6 C , J 7 , and tties (AE) when plotting AE (m^s)- The selection is defined by optimizing 
the signal significance with respect to TZ s i g and TZ^- The solid curve in each plot represents the fit 
projection after correcting for the efficiency of the additional selection (approximately 67% for tttt 
and 88% for Kir). 

The time-dependent CP asymmetries S nn and are determined from a second fit including 
tagging and At information, with the yields and Ak-k fixed to the results of the first fit. The 
At PDF for signal ir + ir~ decays is given by Eq. |l], modified to include the dilution and dilution 
difference for each tagging category, and convolved with the signal resolution function. The At PDF 
for signal Ktt events takes into account B°-B° mixing, depending on the charge of the kaon and 
the flavor of B ta _ g . We parameterize the At distribution in B° — > K + K~ decays as an exponential 
convolved with the resolution function. 

A total of 34 parameters are varied in the fit, including the values of S n7r and C n7T , separate 
background tagging efficiencies for tttt, Ktt, and KK events (12), parameters for the background 
At resolution function (8), and parameters for the background shapes in tties (5), A.E (2), and T 
(5). The signal tagging efficiencies and dilutions are fixed to the values in Table ||, while r and 



Arrid are fixed to their PDG values [11]. To validate the analysis technique, we measure r and Am^ 
in the B Tec sample and find r = (1.66 ± 0.09) ps and Am^ = (0.517 ± 0.062)?i ps _1 . Figure |3| shows 
the asymmetry A mix = (N unmixcd - N mixcd )/ (N unmixcd + N mixcd ) in a sample of events enhanced in 
B — > Ktt decays. The curve shows the expected oscillation given the value of Ara^ measured in 
the full sample. 
The fit yields 

S n7T = -0.01 ± 0.37 (stat)± 0.07 (syst) [-0.66, +0.62], 
C n7T = -0.02 ± 0.29 (stat) ±0.07 (syst) [-0.54, +0.48] . 

For each parameter, we also calculate the 90% confidence level (C.L.) interval taking into account 
the systematic error. The correlation between S wn and C n7T is —13%. Systematic uncertainties 
on S n7T and C n7T are dominated by uncertainty in the shape of the 6 C PDF. Since we measure 
asymmetries near zero, multiplicative systematic errors have also been evaluated (0.05). We sum in 
quadrature multiplicative errors, evaluated at one standard deviation, with the additive systematic 
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Figure 2: Distributions of tbes an d AE (histograms) for events enhanced in signal irir (top) and 
Kir (bottom) decays based on the likelihood ratio selection described in the text. Solid curves 
represent projections of the maximum likelihood fit result after accounting for the efficiency of the 
additional selection, while dashed curves represent qq and irir <-> Kir cross-feed background. 
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Figure 3: The asymmetry -4 m i x between mixed and unmixed events in a sample enhanced in Kir 
decays. The curve indicates the expected oscillation corresponding to Am^ = 0.517 Tips -1 . The 
dilution from qq events is evident in the reduced amplitude near |Ai| = 0. 

uncertainties. Figure |3] shows the At distributions and the asymmetry A wn (At) = (N B o(At) — 
Ng (At))/(N B o(At) + N^ (At)) for tagged events enhanced in signal -kit decays. The selection 
procedure is the same as Fig. |2[ with the likelihoods defined including the PDFs for 8 C , m-Es, 
and AE. 

In summary, we have presented updated measurements of branching fractions and CP-violating 
asymmetries in B° — ► 7r + 7r~, K + ir~, and K + K~ decays. All results are consistent with previous 
measurements. Our measurement of Ak-k is currently the most accurate available, and disfavors 
theoretical models that predict a large asymmetry |L2[ [H| . 
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